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Abstract 

We simulate the improved three-dimensional two-component (ft 4 model 
on the simple cubic lattice in the low and the high temperature phase for 
reduced temperatures down to \T — T c \/T c m 0.0017 on lattices of a size 
up to 350 3 . Our new results for the internal energy and the specific heat 
are combined with the accurate estimates of j3 c and data for the internal 
energy and the specific heat at (3 C recently obtained in ref. pQ. We find 
R a = (1 — A + 1 'A_) / 'a = 4.01(5), where a is the critical exponent of the 
specific heat and A± is the amplitude of the specific heat in the high and the 
low temperature phase, respectively. 



1 Introduction 



The renormalization group (RG) theory of critical phenomena classifies continuous 
phase transitions into so called universality classes, which are characterized by the 
dimension of the system, the range of the interaction and the symmetry of the order 
parameter; see e.g. El EE] At continuous phase transitions, thermodynamic 
quantities follow power laws. E.g. the specific heat behaves as 

C~ A ± \t\- a (l + C±\t\ 9 + ...)+ Cns, (1) 

where a is the critical exponent of the specific heat, 9 the exponent of the leading 
correction to scaling, A + , A_,c + ,C- are amplitudes in the high and low temperature 
phase, respectively. The reduced temperature t = (T — T c )/T c gives the distance 
from the critical temperature T c . c ns is the analytic background, which has to be 
taken into account when a < 0, as it is the case here. Following the RG theory, 
critical exponents are universal, which means that they take exactly the same value 
for any system within a given universality class. Most recent estimates are a = 
—0.0151(3) and 9 = 0.527(13) pQ. For reviews on theoretical and experimental 
results see [Hill]. In addition to the critical exponents, amplitude ratios like A + /A_ 
and c + /c_ are universal, while the individual values of A + ,A_,c + and c_ depend on 
the microscopic details of the model. 

The three-dimensional XY universality class is of particular interest, since the A- 
transition of 4 He is supposed to share this universality class. The experimental study 
of this transition provides by far the most precise experimental results for universal 
quantities like critical exponents and amplitude ratios. Thus this transition gives 
us an unique opportunity to test the ideas of the renormalization group and to 
benchmark theoretical methods. Most recent experiments with 4 He were carried 
out during a spacelab mission Ej- The condition of micro-gravity avoids the 
broadening of the transition due to the gravitational field and hence allows to access 
reduced temperatures down to t ~ 5 x 10 -10 . The most recent analysis [Oj of the 
spacelab data gives a = —0.0127(3) and A + /A^ = 1.053(2) or 

1 - A+/A- , s 

R a = - = 4.154 22 . 2 

a 

Note that R a is much less correlated with the value of a than A + /A_ itself. This 
result can be compared with other experimental results R a = 4.194(19) ^U], R a = 
4.57(40) PU, i? Q = 3.98(2) [12], field theoretic estimates R a = 4.433(77) [13], Monte 
Carlo simulations R a = 4.20(5) Jl] and high temperature expansions combined 
with the equation of state R a = 4.3(2) pQ. Here we make no attempt to give a 
complete overview of theoretical results, we just try to give the most recent and 
hopefully most accurate result for each of the methods. 

Notice that some of these estimates of R a are not consistent among each other. 
This could be interpreted as a violation of universality; however we regard it as more 
likely that systematic errors are underestimated by some of the authors. Here we 
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make an effort to keep errors, in particular the systematic error due to corrections 
to scaling, under control. 

We studied the 4 model on a cubic lattice with periodic boundary conditions 
in each of the directions. The lattice size is L 3 , where L is the linear extension of 
the lattice. The lattice spacing is set to a = 1. The classical Hamiltonian is given 
by 

H=-P'$2?x-$v + Y, + A (<S - X ) 2 ] > ( 3 ) 

<x,y> x 

where the field variable <p x is a vector with two real components. < x, y > denotes 
a pair of nearest neighbour sites. The partition function is given by 



Z((3,X) = J D[0] exp(-H(p,\,<f>)) 



(4) 



where /D[0] denotes the 2L 3 dimensional integral over the field variables. In the 
limit A — > oo the classical XY (plane rotator) model is recovered. Corrections to 
scaling amplitudes such as c± of eq. (JTJ) are functions of the parameter A. It has been 
demonstrated that there exists a value A* at which the amplitudes of the leading 
correction to scaling vanish. The most recent numerical estimate is A* = 2.15(5) 
obtained in ref. pQ. Previous estimates are A* = 2.07(5) and A* = 2.10(6) given in 
refs. [7J [T3] , respectively. 

Here we shall analyse data for the specific heat and the energy density at A = 2.1 
and A = 2.2. In ref. pOj accurate estimates of the inverse critical temperature at 
various values of A are given. In the following we shall use 

P e = 0.5091503(3) [3] at A = 2.1 , 

P c = 0.5083355(3) [4] at A = 2.2 . (5) 

The number in the parentheses gives the statistical error, while the number in 
brackets is an estimate of possible systematic errors. 

The outline of the paper is the following. In the next section we define the 
energy and the specific heat for our lattice model. We summarize the predictions 
of the RG-theory for the free energy density and the specific heat. In section |3] 
we present our numerical results. First we analyse the finite size behaviour of the 
energy density and the specific heat at the critical temperature. Then we discuss 
our results for the low and the high temperature phase. Finally an estimate for R a 
is obtained by fitting these data to the expected power law behaviour. 

2 Critical behaviour of the energy and the spe- 
cific heat 

The so called reduced free energy density is defined by 

f(P,\) = ~]nZ(P,\) , (6) 



2 



which is more convenient for our purposes than the usual — y \nZ. 
We define the energy density as 

E = -% = v T, (7) 

" <xy> 

and the specific heat as 

\ <xy> <xy> 

These definitions differ by factors —1 and (3 2 from standard textbook definitions, 
these factors cancel in the amplitude ratio. Now, let us summarize the predictions 
of the RG-theory for the free energy density. First, the free energy is split into a 
singular and non-singular part 

f = fns + fs ■ (9) 

Let us first discuss the finite size scaling (FSS) behaviour. The singular part of 
the free energy density depends on the parameters of the model and the lattice size 
L. Renormalization group (see e.g. jH3 IT?] ) predicts for a system with periodic 
boundary conditions 

f s (u t , u h , R}, L) = L- d <b(L yt u t , L y »u h , {L^ Ul }) , (10) 

where u t = Ui, Uh = ^2, with i > 3 are the scaling fields (which are analytic 
functions of the Hamiltonian parameters) associated with the reduced temperature 
t {ut ~ t), the magnetic field h {uh ~ h), and the other irrelevant perturbations 
with Hi < 0, respectively, d is the dimension of the system; in our case d = 3. The 
exponent of the thermal scaling field is related with the critical exponent of the 
correlation length: y t = \ jv. Furthermore the hyperscaling relation a = 2 — dv = 
2 — d/yt holds. 

Let us consider the scaling of the energy density and the specific heat at the 
critical point u t = and Uh = 0. Taking the derivative with respect to (3 and 
Taylor-expanding in Ui with i>3we get 

E S (L,(3 C ) = e s L- d+ y>(l + diL ys + d 2 L y * + ...) (11) 

and 

C s {L,p c ) = c s L- d+2 ^{l + fxU 3 + f 2 L* + ...) (12) 

for the singular parts of the energy density and the specific heat, respectively. The 
numerical values of the correction exponents are —ys = uj = 0.785(20) l a and 

-2/4 = uj 2 « 1.8 pap. 
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In the thermodynamic limit, for vanishing external field h = 0, the singular part 
of the free energy density behaves as 

/, = a ± \t\ 2 ~ a (1 + c lilj± |t| A + c li2j± |t| 2A + d x t + c 2il>± |t| A2 + ...) (13) 

with A = vuj and A 2 = vuji- In the following we define the reduced temperature by 
t = (3 — f3 c . Predictions for the singular parts of the energy density and the specific 
heat can be easily derived from eq. (fTHj) by taking the first and second derivative 
with respect to (3. One gets e s = a±\t\ 1 ~ a (1 + ...) and c s = A±\t\~ a (1 + ...) with 

_h — ^± (-iA) 



3 The Numerical results 

Let us first sketch the strategy of our numerical analysis. We parametrise the non- 
singular part of the free energy by its Taylor-expansion in the reduced temperature 
t = P-P e : 

fns(t) = f n8 (0) + e ns t + -C ns t 2 + ... . (15) 

In a first step of the analysis we determine e ns and c ns from the finite size behaviour 
of the energy density and the specific heat at the critical temperature. Then, in the 
analysis of the data for the thermodynamic limit at (3 ^ (3 C , these estimates for e ns 
and c ns are used as input. In this respect, we essentially follow ref. |T3j . 

Note that the amplitudes of the leading correction to scaling are small at A = 2.1 
and A = 2.2. Therefore we shall not take into account the leading correction to 
scaling in our fits. The systematic error introduced this way can be estimated by 
comparing the final results obtained for A = 2.1 and A = 2.2. 



3.1 The energy density and the specific heat at the critical 
temperature 

We analyse data for the energy density and the specific heat obtained on lattices of 
a linear size up to L = 128. The simulations are performed for some (3 S ~ (3 C . We 
have computed the Taylor coefficients of the energy density at (3 S up to 2 nd order. 
This allows us to compute the energy density and the specific heat in a sufficiently 
large neighbourhood of f3 s . 

Most of the data for A = 2.1 were generated already for ref. jTj. Some additional 
data, in particular data for L = 128, were generated more recently for ref. pQ. In 
the case of A = 2.2 the data for L < 16 and for L = 128 were generated for ref. 
PP. Here we have added new data for L = 24,32,48 and 64. In the Monte Carlo 
simulation a mixture of local and wall-cluster ^J] updates was used. For details we 
refer the reader to refs. [3 [T] . 
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First we have analysed the energy density at the central estimates of (3 C given 
in eq. (J3J). In addition, to propagate the error of /3 C , we have analysed the energy 
density at a slightly shifted value for /3 C . We fitted our data using the ansatz 

E = e ns + e s L- d+1 ^ (16) 

without any correction term. We have performed these fits fixing v = 0.6717, which 
is the central estimate of ref. pQ, corresponding to a = —0.0151. In addition we 
performed fits with slightly different values of v to determine the dependence of the 
result for e ns on the input value of v. For v = 0.6717, these fits lead to a x 2 /d.o.f. 
~ 1 starting from L min = 6, where L min is the minimal lattice size that has been 
included into the fit. Furthermore, the result for e ns is very stable when L min is 
increased. Also the statistical error of e ns increases only slowly with increasing 
L m i n . Hence we regard the estimate of e ns obtained from these a fits as reliable. 
Being very conservative, we take our final result from fits with L m i n = 24. Our 
results can be summarized as 

e ns = 0.913213(5) + 20 x ((3 - 0.5091503) + 5 x 10" 7 x (1/a + 1/0.0151) (17) 

for A = 2.1 and 

e ns = 0.913585(5) + 20 x (0 - 0.5083355) + 6 x 10~ 7 x (1/a + 1/0.0151) (18) 

for A = 2.2. For the later use it is convenient to parametrise the dependence of the 
result on the input parameter v in terms of 1/a. 

Next we have fitted the specific heat at f3 c using the ansatze 

C = c ns + c s L- d+2 '» (19) 

and 

C = c ns + c s L- d+2 ^(l + f 2 L-^) (20) 

with v = 0.6717 fixed. Fitting with ansatz (|19p we get x 2 /d-o.f. ~ 1 only for 
Lmin > 16 for both values of A. Hence we performed fits with ansatz (|20|) in 
addition. We have fixed either 0J2 = 1.8 or 002 = 2.0. For this type of fit we get 
X 2 /d.o.f. ~ 1 starting from L min = 6. The results of c ns obtained with 0J2 = 1-8 
and co>2 = 2.0 differ only little. In fig. ^we have plotted the results for c ns obtained 
from the fits discussed above. 

Guided by fig. [TJ we take, for both values of A, our final estimate of c ns from the 
fit with ansatz (|T9|) and L m i n = 24. The error, which is indicated by dashed lines 
in fig. ^ is chosen such that it covers the estimates from fits with ansatz ()19)) as 
well as with ansatz f!20[) . The dependence of the final result on the input values of 
(3 C and v is estimated by redoing the fit with ansatz (|T9"j) and L min = 24 for slightly 
different values of v and (3 C than those chosen above. As result we obtain 

c ns = 157.9(5) + 147000 x (0 - 0.5091503) - 2.1 x (1/a + 1/0.0151) (21) 
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Figure 1: The non-singular part of the specific heat c ns obtained from numerical 
data at (3 C as a function of the smallest lattice size L m i n that has been included in 
the fit. We have used either an ansatz without any correction to scaling (jl9j) or an 
ansatz with sub leading corrections to scaling ()20|) with an RG-exponent 2/4 = —1.8. 
For details see the text. 
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for A = 2.1 and 

c ns = 155.6(4) + 121000 X (/3 - 0.5083355) - 2.1 x (1/a + 1/0.0151) (22) 
for A = 2.2. 

3.2 The energy density in the thermodynamic limit 

The simulations in the high temperature phase were already discussed in ref. pQ. In 
the high temperature phase we expect that the observables converge exponentially 
fast to the thermodynamic limit as L — > oo. Throughout we have used L > 10£ in 
these simulations. In particular for our largest correlation length £ pa 30 we used 
L = 350. Hence within our numerical accuracy, the results should coincide with the 
thermodynamic limit. 

In the low temperature phase, the breaking of the continuous U(l) symmetry 
leads to a Goldstone mode. As a result, leading corrections to the thermodynamic 
limit of the energy density are 0(L~ 3 ) [SHI |2l]. We fitted our data for the energy 
density with the ansatz 

E(L) = E(oo)+cL- 3 . (23) 

Typically, the difference E(oo) — E(L max ) of the fit result for the thermodynamic 
limit and the result for the largest lattice size L max that we have simulated is of 
a similar size as the statistical error of E(L max ) and E(oo). Therefore we are 
confident that the quoted error, which is the statistical error of the result of the fit 
with ansatz (|23)1 . is reliable. For illustration we have plotted in Fig. (J2J) the Monte 
Carlo results for the energy density E at A = 2.1 and (3 = 0.51 as a function of L~ 3 . 
Our final results for the thermodynamic limit of the energy density are summarized 
in table H in addition we give the second moment correlation length ^2nd in the 
high temperature phase. 

We have also measured the specific heat. However the result for R a obtained 
from these data is consistent with but less precise than that obtained from the 
energy density. Therefore we skip the discussion of the specific heat data. 

First we have fitted our data for the energy density with the ansatz 

E = e ns + c ns {(3 - p c ) + a±\/3 - (3 c \ l ~ a . (24) 

In these fits, we take a, (3 C and the corresponding values of e ns , eqs. ()17I18J) . and 
c ns , eqs. fl2H22|) . as input. 

It turns out that in the case of A = 2.1, fits that include \/3 — f3 c \ < 0.005 have 
a large x 2 /d.o.f. . To understand this problem we resorted to a simpler analysis of 
the data. We took pairs of inverse temperatures such that (3 C — j3\ pa /? 2 — /?c - These 
two values are sufficient to determine the parameters a + and a_ of the ansatz (|24|). 
In Fig. we give the results for R a obtained this way as a function of (/% — /?i)/2. 
We see that the estimate of R a is roughly linear in f3 2 — Pi- A linear extrapolation 
to p 2 — Pi — suggests a value of R a that is slightly larger than 4. 
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Figure 2: The energy density E at A = 2.1 and (3 = 0.51 is plotted as a function of 
L~ 3 , where L is the linear lattice size. The circles give our Monte Carlo results for 
the lattice sizes L = 128, 144, 192, 256 and 288. The triangle gives the result of the 
extrapolation to L = oo with ansatz (|23p. In the fit L — 128 had been excluded. 
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Table 1: Results for the energy density E of the 4 model at A = 2.1 and A = 2.2 
in the thermodynamic limit for various values of the inverse temperature (3. In 
addition we give second moment correlation length £ 2n d m the high temperature 
phase. 
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Figure 3: Estimate of R a obtained from two values of the inverse temperature (3\ 
and /?2 using the ansatz (ffifj) . Only results for A = 2.1 are given. Two errors are 
displayed: The smaller one is the due to the statistical error of the energy density 
at /?i and f3 2 ] the larger one is due to the uncertainty of c ns , e ns and [3 C . The dotted 
line should only guide the eye. 
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Motivated by this observation, we performed fits with the ansatz 



E 



e ns + c ns (P - p e ) +a ± \P- (3 c \ l ~ a + d{(3 - (3 C ) 2 , 



(25) 



where we have added a further term to the Taylor expansion of the analytic part 
of the energy density. Note that the first analytic correction to the singular part 
comes with a very similar exponent: 2 — a. Hence the fit parameter d will provide 
only an effective amplitude for the combination of the two terms. Note that also 
the exponent of subleading corrections 1 — a + A2 is only slightly larger than 2. We 
also tried to explicitely take into account these terms in the fit. This leads however 
to very large errors for the coefficients. 

Using ansatz (|23jl we get fits with a x 2 /d.o.f. smaller than one for the interval 
0.506 < /3 < 0.512 and 0.5055 < p < 0.511 for A = 2.1 and A = 2.2, respectively. 
Using the central values of c ns , e ns , p c and a as input we get R a = 4.011(5) and 
R a = 4.017(6) for A = 2.1 and A = 2.2, respectively. Remember that for A* = 
2.15(5) leading corrections to scaling vanish [1. Since the results for R a at A = 2.1 
and A = 2.2 are almost the same, we conclude that leading corrections to scaling 
can be safely ignored. 

In order to check for the effect of corrections to the ansatz (|23|) discussed above, 
we have repeated the fits using /3-values that are further off from P c . In particular, 
fitting the data for p = 0.503, 0.505, 0.506, 0.512, 0.513, 0.515 in the case of A = 2.1 
we obtain R a = 4.006(3) and fitting the data for p = 0.501,0.5035,0.5055,0.511, 
0.512 in the case of A = 2.2 we obtain R a = 3.988(2). Note that the average of 
\P — p c \ for this second set of fits is more than twice as large as for the first set of 



Next we have computed the error of R a due to the uncertainty of the input 
values of c ns , e ns and P c . To this end, we have repeated the fits using the central 
values of c ns , e ns and P c shifted by the error estimate. Summing the errors of R a 
due these input values we get a little less than 0.02 for both values of A. 

Finally, we have also repeated the fits for a 7^ 0.0151 to obtain the dependence 
of our numerical estimate of R a on a. 

Summing up all errors discussed above, we arrive at the final estimate 



Notice that the dependence on a is rather small; e.g. inserting the estimate a = 
-0.0127 of ref. we get R a = 3.99. 

4 Conclusions and Comparison with the litera- 
ture 

We have studied the universal amplitude ratio A+/A- of the specific heat in the 
three dimensional XY universality class. Since |a| is rather small, it is difficult 



fits. 



R a = 4.01(5) - 8 x (a + 0.0151) 



(26) 
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to disentangle the singular and the non-singular part of the specific heat or the 
energy density. This problem holds for numerical data obtained from Monte Carlo 
simulations of lattice models as well as for experimental data. 

While in ref. jHj the standard XY model was simulated we have studied the im- 
proved (p 4 model, allowing us to ignore leading corrections to scaling in the analysis 
of the data. Given the problem discussed above, this is an important advance. 

In our analysis we have combined information from the finite size scaling be- 
haviour at the transition [7j with precise results for the thermodynamic limit 
in the low and the high temperature phase. Reaching correlation lengths up to 
£,2nd ~ 30 subleading corrections to scaling are under control. 

We made an effort to take carefully into account various sources of systematic 
error. Pessimistically we have summed up these errors to arrive at our final estimate 
R a = (1 - A + /A_)/a = 4.01(5). 

Our estimate for R a is significantly smaller than most theoretical and experimen- 
tal results. There is only good agreement with the experimental result R a = 3.98(2) 
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